Robustness of interdependent networks under targeted attack 
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Abstract 

When an initial failure of nodes occurs in interdependent networks, a cascade of failure between 
the networks occurs. Earlier studies focused on random initial failures. Here we study the robust- 
ness of interdependent networks under targeted attack on high or low degree nodes. We introduce 
a general technique and show that the targeted- attack problem in interdependent networks can be 
mapped to the random-attack problem in a transformed pair of interdependent networks. We find 
that when the highly connected nodes are protected and have lower probability to fail, in contrast 
to single scale free (SF) networks where the percolation threshold p c = 0, coupled SF networks 
are significantly more vulnerable with p c significantly larger than zero. The result implies that 
interdependent networks are difficult to defend by strategies such as protecting the high degree 
nodes that have been found useful to significantly improve robustness of single networks. 
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Modern systems due to technological progress are becoming more and more mutually 
coupled and depend on each other to provide proper functionality pj-yj. Social disruptions 
caused by recent disasters, ranging from hurricanes to large-scale power outages and terrorist 
attacks, have shown that the most dangerous vulnerability is hiding in the many interdepen- 
dencies across different networks 4 L The question of robustness of interdependent networks 
has recently become of interest {5 -a]. In interdependent networks, nodes from one network 
depend on nodes from another network and vice versa. Consequently, when nodes from one 
network fail they cause nodes in the other network to fail too. When some initial failure of 
nodes happens, this may trigger a recursive process of cascading failures that can completely 
fragment both networks. 

PI 

Recently, a theoretical framework was developed [7] to study the process of cascading 
failures in interdependent network caused by random initial failure of nodes. They show that 
due to the coupling between networks, interdependent networks are extremely vulnerable 
to random failure. However, when we consider real scenarios, initial failure is mostly not 
random. It may be due to a targeted attack on important central nodes. It can also occur to 
low central nodes because important central nodes are purposely defended, e.g. in internet 
networks, heavily connected hubs are purposely more secured. Indeed, it was shown that 
targeted attacks on high degree nodes 9l4l3|] or high betweeness nodes [14j in single networks 
have dramatic effect on their robustness. The question of robustness of interdependent 
networks under targeted attack or defense has not been addressed. 

In this Letter, we develop a mathematical framework for understanding the robustness 
of interdependent networks under initial targeted attack which depends on degree of nodes. 
The framework is based on a general technique we develop to solve targeted attack problems 
in networks by mapping them to random attack problems. A value W a (ki) is assigned to 
each node, which represents the probability that a node % with ki links is initially attacked 
and become inactive. We focus on the family of functions {12I ] 

k a 

W a (ki) = * -00 < a < +00. (1) 

When a > 0, nodes with higher degree are more vulnerable and those nodes are intention- 
ally attacked, while for a < 0, nodes with higher degree are defended and so have lower 
probability to fail. The case a = 0, Wq = -k, represents the random removal of nodes [t| 
and the case a — > 00 represents the targeted attack case where nodes are removed strictly 
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in the order from high degree to low degree. For the a < case, nodes with zero degree 
should be removed before analysis begins. 

Our model consists of two networks, A and B, with the same number of nodes, N. The 
N nodes in each network are connected to nodes in the other network by bidirectional 
dependency links, thereby establishing a one-to-one correspondence. The functioning of a 
A-node in network A depends on the functioning of the corresponding B-node in network 
B and vice versa. Within each network, the nodes are randomly connected with degree 
distributions Pa(^) and Pb(/c) respectively. We begin by studying the situation where only 
network A is attacked. We initially remove a fraction, 1 — p, of the A-nodes of network A 
with probability W a (ki) (Eq.flTJ) and remove all the A-links that connect to those removed 
nodes. As nodes and links are sequentially removed, network A begins to fragment into 
connected components. Nodes that are not connected to the giant component are considered 
inactive and are removed. Owing to the dependence between the networks, all the B-nodes 
in network B that are connected to the removed A-nodes in network A are then also removed. 
Network B also begins to fragment into connected components and only the nodes in the 
giant component are kept. Then network B spreads damage back to network A. The damage 
is spreaded between network A and B, back and forth until they completely fragment or 
arrive to a mutually connected component and no further removal of nodes and links occurs. 

The main idea of our approach is to find an equivalent network A', such that the targeted 
attack problem on interdependent networks A and B can be solved as a random attack prob- 
lem on interdependent networks A' and B. We start by finding the new degree distribution 
of network A after removing, according to Eq.flTJ, l—p fraction of nodes but before the links 
of the remaining nodes which connect to the removed nodes are removed. Let A p (k) be the 
number of nodes with degree k and P p (k) be the new degree distribution of the remaining 
fraction p of nodes in network A, 

m) = (2) 

When another node is removed, A p (k) changes as 

A ip -y N) {k) = A p {k)- (3) 
where < k(p) a >= ^P p (k)k a . In the limit of N — > oo, Eq.flS]) can be presented in terms of 
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derivative of A p (k) with respect to p, 



dA p (k) P p (k)k a , . 

dp < k(p) a >' 

Differentiating Eq.(j2]) with respect to p and using Eq. (J5J), we obtain 

dP p (k) _ P P (k)k a 
V ^p~ ~ FpW <ife(p)«>' { ) 

which is exact for N — > oo. In order to solve Eq.([5]), we define a function G a (x) = 

^2 h P(k)x ka , and substitue / = G~ 1 (p). We find by direct differentiation that 

satisfy the Eq.([5]). With this degree distribution, the generating function of the nodes left 
in network A before removing the links to the removed nodes is 

G Ab (x) = P P (k)x k = -£ P(k)f a x k . (8) 

k P k 

On the other hand, the fraction of the orginal links that connect to the nodes left is 

_ = pN < k{p) > E k P(k)kf ka 

P ~ N<k> E k P{k)k ' 1 ' 

where < k > is the average degree of the original network A, < k(p) > is the average degree 
of remaining nodes before the links that are disconnected are removed. So the generating 



function of the new degree distribution o: 



the removed nodes are also removed is 21] 



the nodes left in network A after their links to 



G Ac (x) = G Ab (l~p + px). (10) 

The only difference in the cascading process under targeted attack from the case under 
random attack is the first stage where the initial attack is exerted on the network A. If we 
find a network A' with generating function Gao(x), such that after a random attack with 
(1 — p) fraction of removed, the generating function of nodes left in A' is the same as GU c (x), 
then the targeted attack problem on interdependent networks A and B can be solved as a 
random attack problem on interdependent networks A' and B. We find Gao{x) by solving 
the equation Gao(1 — p + px) = Ga c {x) and from Eq. lflQj) . 



G A0 {x)=G Ab (l + P -(x-l)). (11) 



Up to now, we have mapped the problem of cascade of failures of nodes in interdependent 
networks caused by initial targeted attack to the problem of random attack. We can see that 
the evolution of equations only depends on the generating function of network A, but not on 
any information about how the two networks interact with each other. Thus, this approach 
can be generally applied to study both single networks and other interdependent network 
models. 

Next we can apply the framework developed in Ref. Q], qa{p) = 1 — <5ao[1 — p(l — /a)], 
where Ja is a function of p that satisfies the transcendental equation Ja = GUi[l — p(l — /a)]- 
Analogous equations exist for network B. As the interdependent networks achieve a mutually 
connected giant component, the fraction of nodes left in giant component is p^. The system 
satisfies the equations 

x = pg A (y), 
V = P9b(x), 

where the two unknown variables x and y satisfy p^ = xgg(x) = ygA^y)- Eliminating y 
from these equations, we obtain a single equation 

x = P9a[p9b(x)}. (13) 

The critical case (p = p c ) emerges when both sides of this equation have equal derivatives, 

which, together with Eq.( TT3|) . yields the solution for p c and the critical size of the giant mu- 
tually connected component, Poo(Pc) = x c gB(x c ). In general, there is no explicit expression 
as a solution and p c and x c can be found numerically. 

We now analyze the specific classes of Erdos-Renyi (ER) 15|, ll6[ and scale-free (SF) 



19| networks. The lines in Fig. [T] represent the critical thresholds, p c , for coupled coupled 
SF networks with different a obtained by solutions of Eq. (1131) and Eq. (fl4j) . which are in 
excellent agreement with simulations. Several conclusions from Fig. [1] are as follows: (i) 
Remarkably, while p c for a single SF network approaches to quickly when a becomes 
zero or negative (see also fl2j|). p c for interdependent networks is non-zero for the entire 
range of a (Fig. (Ha)). This follows from the fact that failure of the least connected nodes 
in one network may lead to failure of well connected nodes in the other network, which 
makes interdependent networks significantly more difficult to protect compared to a single 



network, (ii) targeted attacks (a > 0) and defense strategies (a < 0) are more effective 
for interdependent networks with broader degree distributions. In Fig. [T](b), comparing the 
lines of A = 2.5, A = 2.8 and A = 3.4 with m = 2, one can see that the lower is A the more 
sensitive is p c to the change of a. Accordingly, robustness of interdependent networks with 
broader degree distributions decreases more under the same targeted attacks. 

Simplified forms for G A b(x) , G Ac (x) and G A q{x) from Eqs.(!8l. (jT0l) and (TTTT) exist when 
a = 1, 

G Ab {x) = ~J2 P(k)f k x k = -G A0 (fx), (15) 

v T v 

G Ac {x) = l -G M {f{l-p + px)) 1 (16) 

G A0 (x) = -G A0 ^f(x-l) + f). (17) 
p p 

where G A q(x) is the original generating function of the network A, / = G~j^{p) and p = 

Explicit solutions of percolation quantities exist for the case of interdependent Erdos- 
Renyi networks, when a = 1 and both of the two networks are initially attacked simuta- 
neously. The two networks originally have generating functions G A q(x) and Gbo{x). Initially, 
(1 — pi) and (1 — p 2 ) fraction of nodes are targeted (according to Eq. (1) and a = 1) and 
removed from network A and B respectively. Similarly, we start by finding the equivalent 
networks A' and B' such that a fraction (1 — p\P2) of random initial attack on both of the 
networks has the same effect as (1 — px) and (1 — p%) fraction of nodes are intentionally 
removed from network A and network B respectively. After removal of initially failed nodes 
and all the links that connect to the removed nodes, according to Eq.( fl6l) . the generating 
function of the nodes left in network A is 

G Ac (x) = -^G A0 (fx(l-px+Pix)), (18) 

1 G' ( f ) 

where fx = G^Jpi), p\ = /i A , ■ Furthermore, (1— p-z) fraction of the remaining A-nodes 
are randomly removed. Because each remaining A-node's corresponding B-node in network 
B has a possibility (1 — pi) to be initially attacked, which leads to fail this A- node. The 
generating function of the nodes left in network A is 

G Ad {x) = G Ac (l -p 2 + P2x) = —G AQ (fi + pxfxP%{x - 1))- (19) 

Pi 



6 



Now we can find the generating function of the equivalent network A' by G A q(1 — P1P2 + 
Pip 2 x) = G Ad (x): 

G A0 (x) = -G A0 (^h(x - 1) + h). (20) 
Pi Pi 

The same holds for network B' . 

For ER networks, the generating function is G (x) = e < k >( x - 1 ) j^lj^ so f x = + 1,^2 = 
§g + 1, G Ao (x) = G A1 (x) = e <k> ^ x ~V and G B0 {x) = G m {x) = e <k> ^ x ~ l \ From 
Eq.tfTJ, 

X = PlP29A{y) = PlP2(l - /a), ^ 

V = PiP29b(x) = Pip 2 (l - /b), 



where 

j A = e <fc>i/ 1 2 ^(/ j4 -i) ) 
f B = e <fc>a^*(/B-i) > 
In the case < >i=< >2=< > and pi = p 2 = P, we find that 



(22) 



Poo =p 2 (l-e <k>f2p °°). (23) 



where fi = f 2 = f = + 1, and p c satisfies relation: 



<fc> 

,2 



< k > p 2 J c = 2.4554, (24) 

with f c = lr ^ly + 1. Fig. |2] shows that the simulation confirms well the theory. Compared 
to the case of random attack on one network, where p c = 2.4554/ < k > [7j, in Eq.f l24|) . 
the factor f c reflects the effect of targeted attack on high degree nodes to increase p c . 
The term p 2 c in Eq. (l24"]) is since we are initially attacking both networks simutaneously 
instead of only attacking one network. Indeed for the case of initial random attack on 
two networks simultaneously, from Eq.( l2T|) and f A = e <k>iy ^ A ~ 1 \ fs = e <k>iy ^ B ~^ [7| we 
obtain < k > p 2 c = 2.4554. 

In summary, we developed a theoretical framework for understanding the robusteness of 
interdependent networks under targeted attacks on specific degree nodes. We introduce a 
method and show that targeted-attack problems in networks can be mapped to random- 
attack problems by transforming the networks which are under initial attack. It provides a 
routine method (if the random-attack case is solvable) to study the targeted-attack problems 
in both single networks and randomly connected and uncorrelated interdependent networks, 
i.e. (i) the case of three or more interdependent networks, (ii) the case of partially coupled 



interdependent networks, (iii) the case in which a node from network A can depend on more 
than one node from network B. By applying the method, we find that in contrast to single 
networks, when the highly connected nodes are defended (a < 0), the percolation threshold 
p c has a finite non-zero value which is significantly larger than zero. For example, when 
the degrees of all nodes are known and nodes can only be damaged from lower degree to 
high degree (a — > — oo), p c ~ 0.46 for coupled SF networks with A = 2.8 and < k >= 4 
while p c for the same single SF network is (Fig. [1]). The implications of the present study 
are dramatic. The current methods applied to design robust networks and improve the 
robustness of current networks, i.e. protecting the high degree nodes, need to be modified 
to apply to interdependent network systems. 
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FIG. 1: (a) Dependence of p c on a for SF single and interdependent networks with average 
degree < k >= 4. The lower cut-off of the degree is m = 2. The horizontal lines represent 
the upper and lower limits of p c . The black dashed line represents p c for SF free network, (b) 
Values of p c vs a for SF interdependent networks with different A and lower cut-off m = 2. 
The A in the legends of both the graphs are approximate numbers. 



10 




p 

FIG. 2: Values of p,^ vs p when both networks are initially attacked. Both networks in the 
interdependent networks are ER or SF networks with the same average degree. The symbols 
represent simulation data (N = 10 6 nodes). The solid lines are theoretical predictions, 
Eq.( l23p . The dashed line represents simulation data for interdependent scale-free networks 
with A = 2.8, < k >= 4. All results are for a = 1. Inset: Values of p c vs average degree of 
ER networks. The symbols represent simulation data, while the line is the theory, Eq. (|24p . 
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